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VQ We use the threshold expansion and non-relativistic effective theory to determine 

^^_ ^ the bottom quark mass from moments of the bb production cross section at next- 

^^ , to-next-to-leading order in the (resummed) perturbative expansion, and including 

lyL I a summation of logarithms. For the MS mass Wit, we find rni,(m(,) = (4.26 it 

^ ■ 0.12) GeV. 

a^ 

--^ ■ 1 Introduction 

^ i' In this talk we report on a determination of the bottom quark mass at next-to- 

next-to-leading order (NNLCX) based on dispersion relations for e+e^ -^ bbX. 
This idea dates back to Ref. El and relies on the relation 
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n^(l)U^O= 31+1 ^''''(^)' (1) 



between the two-point function of the bottom vector current and the inclusive 
5^ I bottom quark cross section mediated by a virtual photon, which is valid up to 

a small correction from bb radiation from light quarks. The left hand side can 
be computed in perturbation theory; the right hand side from data. 

The parameters of the lowest T{nS) resonances are well-measured, but the 
continuum cross section above the threshold is not. Hence the experimental er- 
ror of the right hand side decreases with increasing n, because higher moments 
weight lower s. When the integral over s is saturated by the threshold region, 
the bottom quarks are non-relativistic, and 2m^,l ^/n and nib/n appear as new 
momentum scales in the problem, related to the typical momentum mf,v and 

"^Talk presented at RADCOR98, Barcelona, September, 1998. The numerical result of this 
version supersedes the version printed in the proceedings volume, which was affected by a 
program error. 



non-rclativistic energy mi,v'^ of the bottom quarks, respectively. When n is 
large enough such that v ^ as, ordinary perturbation theory breaks down, be- 
cause there are terms of the form {ag/v)^ to all orders in perturbation theory. 
Even when v is only small, but not as small as a^, one may want to resum 
those enhanced terms systematically. 

Defining w = (1 — 4m^/s)^/^, where rri}, is the h quark pole mass, the cross 
section e+e^ ^ bbX for s « 4?Ti^ can be expanded in a double series in ag and 
v. 

Rbb^vY^ Y^ Cki a'; u' X logs of V. (2) 

k=0 l=-k 

In the resummed perturbative expansion, a N^LO calculation has to account 
for all terms with k + I < p. Such a rearranged expansion has Jirst been 
considered in this context in Ref. □ to LO accuracy and in Ref. u to NLO 
accuracy. The step to NNLO accuracy implies new difficulties, because at this 
order the naive relativistic or non-relativistic approximations lead separately 
to divergent integrals. One has to be more precise about factorizing the two 
momentum regimes. 

Powers of v in (0) arise from ratios of momentum scales. We have to dis- 
entangle the contributions from the different scales in order to be sure that for 
a high-order loop graph, which cannot be calculated exactly, we have taken 
into account all terms with k + I < 2. In Ref. u we described how to expand 
a Feynman integral in v to any given order, without calculating the exact 
result. The expansion procedure is based on decomposing the loop integral 
into contributions from different loop momentum regions and has a transpar- 
ent interpretation in terms of effective theories. The present application, the 
perturbative calculation of i?;,;, in the threshold region, has much in common 
with a QED bound state calculation done with non-relativistic QEDEfQ. The 
main difference is that we perform factorization in dimensional regularization, 
rather than using momentum space cut-offs and a photon mass as infrared 
regulator, and take advantage of the fact that scaleless integrals vanish in di- 
mensional regularization. The absence of any explicit cut-off scale leads to 
a homogeneous expansion in v, manifest power counting and makes analytic 
calculations much easier at two loops. The unfamiliar dimensional renormal- 
ization of the Schrodinger problem is a negligible price to pay in comparison 
to these advantages, because it is not necessary to obtain the Green function 
exactly in d dimensions. 

In the foUowing^e summarize the calculation of e+e^ — > bbX at NNLO 
using these methodsLla and we determine the bottom quark mass by calculating 
the moment integral (M. In addition to the NNLO terms in the sense of (0), 
we also sum all logarithms of v of the form (aslnw)'^ and ag {aglnv)''. A 



subset of NNLO logarithms of the form a^ (aslnv)'' can easily be summed 
as well, but a complete summ.ation requires another difficult calculation to 
be done. This write-up concentrates on necessary formulae. A more general 
perspective on the problem, and Ji, more complete list of references related to 
the subject can be found in E^LrLWhile we were preparing this work, other 
NNLO calculations appearedE30'E3. A brief comparison with these results is 
given at the end. 

2 Heavy quark production cross section near threshold 

2. 1 Classification of momentum regions 

The expansion (Q) can be constructed by decomposing each Feynman integral 
that contributes to the vacuuxp polarization n(q^) into contributions from the 
following momentum regions cl: 



hard (h) 

soft (s) 

potential (p) 

ultrasoft (us) 



/o ^ riib, I ^ mb, 

Iq ^ mbV, I ^ mi,v, 

lo ^ rnbV^, I ^ nibV, (3) 

Iq ^ iTibv'^, I ^ nibv'^, 



assuming a frame where q = 0. In dimensional regularization, the factorization 
of the different momentum regions is achieved by appropriate expansions of 
the Feynman integrands, without the need for explicit cut-offs, so that one can 
integrate over the whole integration domain in each region. As a consequence 
of this expansion, all terms have a definite homogeneous scaling behaviour in 

V. 

In the following, we integrate out first the hard, relativistic modes. This 
defines the NRQCD Lagrangian in dimensional regularization. In a second 
step, we also integrate out soft modes and potential gluon modes to arrive at 
another effective theory, 'potential NRQCD' (PNRQCD), the QCD analogue 
of PNRQED, introduced in Ref. t3. The final result is obtained by computing 
the vacuum polarization in PNRQCD perturbation theory, with the Coulomb 
interaction treated non-perturbatively 

2.2 Relativistic corrections 

At tree-level relativistic corrections arise from the anti-particle pole of the 
heavy quark propagator and the anti-particle components of the heavy quark 
four-spinor field. Beyond tree level, the hard loop momentum regions induce 
further terms in the non-relativistic effective theory. 



The heavy quark current correlation function n(g^) contains IPI, hard 
subgraphs which (a) do not connect to the virtual photon vertex, (b) connect 
to one of the virtual photon vertices and (c) connect to both virtual photon 
vertices. The last possibility is of no interest, since it cannot contribute to the 
discontinuity of Il{q^) and hence is irrelevant for the computation of {^. 

Acrniintipg for the hard subgraphs (a) leads to the familiar NRQCD La- 
grangianliJ^E£l. At NNLO, the following terms are needed: 

/:nrqcd = ^t hi?o + ^ j V- + g^ ^t^V - ^ V-V • Bi^ 

+ -^ i^Uo- E- E- D)tlj + -^ i;^a-(DxE-ExD)tP 
8m^ V / 8m^ V / 

+ antiquark terms +£iight- (4) 

Since we work in dimensional regularization, etjk should be defined carefully. 
As long as we are only interested in n(q^), we can avoid the problem by 
writing the Lagrangian in terms of anti-commutators of Pauli-matrices. Hence 
a ■ B must be interpreted as i/4 [a', cr^] G*-'. A similar interpretation holds 
for the spin-orbit interaction. The scaling of the interaction terms relative to 
ipHD^ip depends on whether the quark and gluons fields are considered as soft, 
potential or (in the case of gluons) ultrasoft. For potential quarks, D'^/{2mb) is 
a leading order term, so its coefficient would be needed to order a^. However, 
since the coefficient is 1 to all orders in PT, there is nothing to calculate. The 
chromomagnetic interaction is suppressed by only one power of v for a soft 
giuon and by two powers of v for an ultrasoft gluon. However, since one needs 
at least two chromomagnetic insertions to obtain a non- vanishing contribution 
to H(g^), it is sufficient to know di at tree level. Since, for the time being, we 
sum logarithms of v only at NLO and not at NNLO, we have di ~ d2 ^ d^^ ~ 1. 
Accounting for the hard subgraphs (b) leads to the effective non-relativistic 
j*bb coupling 

Q7*Q = ci V' V^x - A V' V^ (iDfx + ■■■, (5) 

where the ellipsis refers to terms not needed for H(g^) and at NNLO. At NNLO, 
we can use C2 = 1, while ci is needed at order a^. The required matching 
calculation has been done in Refs.Ej'tJ. The result is 
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where 



and 60 = 11 ^ 2n^/3. The term proportional to Si sums all next-to- leading 
logarithms of the form as{as Inw)''. There are no leading logarithms and this 
term is the only source of next-to-leading logarithms in the problem. 

2.3 Instantaneous interactions 

The loop diagrams constructed from the NRQCD Lagrangian contain soft, 
potential and ultrasoft modes. We, now integrate out the soft modes and 
potential gluons. As discussed inQil3, these modes give rise to instantaneous 
interactions. At NNLO, it is sufficient to match the (on-shell) quark-antiquark 
scattering amplitude to the required order. We do this order by order in as, 
and, except for 02 below, we have checked that the resulting potentials are the 
same whether we use Coulomb gauge or a general covariant gauge. The terms 
in the effective PNRQCD Lagrangian, which we need at NNLO, are given by 



2mf, 8ml j \ 2mb &ml 



/:PNRQCD = V^M *9° + T^ + ^ 1 ^ + XM *50 
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where we indicated the spin and colour indices. In general, the PNRQCD La- 
grangian contains local interactions of potential quarks and ultrasoft gluons, 
which give rise to retardation effects, but they contribute only at N'^LO, pro- 
vided one counts asimtv^) as v, as one does for asimtv) and asinib). Power 
counting confirms the well known fact that for r ^ l/{mbv) and as ^ v, 
the leading order Coulomb interaction is not suppressed relative to the free 
heavy quark Lagrangian. The instantaneous interaction in the third line can 
be treated perturbatively. At NNLO, we need all potentials of the form a^'^/r, 
a^/r"^ and ag/r^, counting 5^^''{r) as 1/r^. To compute Il{q^), only the colour 
singlet and spin-1 projection is needed. The result, obtained from matching 
the quark-antiquark scattering amplitude, is, in momentum space, 

SVijkl,abcd{p,q) = SijSkl6abScd ■ — r^ • 
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r(i 



where 6i = 102 — 38n//3 and ai and 02 denote the one-loop and two-loopE3 
radiative corrections to the Coulomb potential, respectively. The potentials 
more singular than 1/r lead to ultraviolet divergent integrals in PNRQCD 
perturbation theory. These divergences cancel with divergences that arise in 
the calculation of ci at the two-loop order. To obtain the correct finite terms, 
it is necessary to compute the potentials that lead to divergent insertions to 
order e, where d = 4 — 2e is the number of space-time dimensions. 



2.4 Perturbation theory with potential NRQCD 

The heavy quark current correlation function is now obtained from the two- 
point functions of the effective currents <M), computed with the PNRQCD 
Lagrangian. Because the unperturbed PNRQCD Lagrangian contains the 
Coulomb interaction, the propagator for a bb pair is the Coulomb Green func- 
tion Gc{x,y;E), where E = \fq^ — Irrib. At LO in PNRQCD perturbation 
theory, one finds 



n(g^ 



2ml 



Gc(0,0;S). 



(10) 



The Creen function at the origin is ultraviolet divergent. We need its value 
in dimensional regularization, with MS subtractions. To this end we write all 
integrals first in monrcntum space and note that all divergences can be removed 
by counterterms that involve only a finite number of loops. In the case of 
60(0, 0; i?), for example, all diagrams with two or more giuons exchanged are 
finite. The result in the MS scheme is 



Gc(0,0;£;) = 



■in 



1 1 , -Am^E 1 
^+2ln— ^--+,.+^(l-A) 



(11) 



where A = 2as / {'i\J ~ E / mb) . The constant —1/2 is scheme-dependent. Since 
the cross section requires only the discontinuity of H, the subtraction procedure 
is irrelevant for the moments. However, essentially the same procedure can be 
used to obtain more complicated integrals in dimensional regularization, which 
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are needed at NNLO. In particular, the pole parts of the divergent potential 
insertions are all proportional to Gc(0, 0; E) as they have to be in order that 
the /i dependence that comes from (g) cancels at NNLO. 

To obtain the final result at NNLO one has to compute: (i) a single in- 
sertion of the kinetic energy correction 9^/(8to^) to the quark propagator; (b) 
the single and double insertion of the terms in dV{p, q), which are suppressed 
by one power of w ^ a^ compared to the LO Coulomb potential; (c) the sin- 
gle insertion of all terms in 5V{p,q), which are suppressed by two powers of 
V ^ Us- The result is an expression for R^^i that can be cast into the form 

i?,S = a,{/o(A,0+as/i(A,0 + «?/2(A,0}, (12) 

where I — ln(— 4?n(,i?//x'^). One can then compute the integral (|l|) and compare 
it with its experimental value. Note that R^i contains a continuum starting at 
4?7i^ and an infinite series of bb bound state poles, which are included in the 
integral. In (^2|) the denominators of the bound state pole contribution are 
expanded around the LO pole position. Before we turn to a numerical analysis 
of the moments of (|l2|), we have to discuss the issue of mass rcnormalization. 
We also mention the following checks we performed on ([I4). We re- 
expanded (|l2| ) to order a1 and found agreement with the cross section near 
threshold computed to this order in Ref.EZl. We also obtain the bound S|tate pole 
position and its residue in analytic form and confirm the result of Ref.ll3. Both 
together are strong checks that the factorization in dimensional regularization 
has been done correctly. Cojiversely, our calculation provides an independent 
check of Lhe_Eesult of Ref. tZi, which has been used in the other NNLO cal- 
culations t2Hlilll3, which followed the 'direct matching' procedure suggested in 
Ref.llj, rather than factorization in dimensional regularization. 

3 Potential subtracted quark mass 

The heavy quark production cross section near threshold is conventionally 
expressed in terms of the quark pole mass. If one uses another mass rcnor- 
malization convention that differs from the pole mass by 5m, one finds terms 
of the form (Sm/E)'', which modify the structure of (||), and whirb sppm to 
complicate the resummation. However, the pole mass is known Ej'lJ to be 
more infrared sensitive than the heavy quark production cross section itself 
and one therefore expects that a badly convergent series expansion for Rf^j^, 
when expressed in terms of mi,, would prevent us from extracting mi, accu- 
rately. The convergence should be improved, when R/^i is expressed in terms 
of a less IR sensitive mass parameter. To implement this observation in the 
calculation, we make use of a systematic cancellation of infrared contributions 



to the pole mass and the Coulomb potential in coordinate-spaceE^iiJ and define 
the potential-subtracted (PS) quark mass mt,^ps{l^f) byE3 
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Explicit expressions for Srni^2 can be found in Ref.E3. Note that m;, — ?7if,,ps(/^/) 
is proportional to a subtraction scale /i/, which must be chosen to be smaller 



than rubV. We insert (13) into (h2|) and expand the small correction terms 
involving 6mi^2- However, the term A = 4as/i//(37r) is not expanded, when 
rrib is replaced in E, X or I, because A, which counts as order mbV^ , is of the 
same order as E = ^/q^ — 2mb. The result is an expression of the same form 
as (p^), but with mb^psil^f) ^^ input parameter. If our expectation is correct, 
the expansion (n% should be more convergent in this new variable. We then 
extract mb^psip-f) from comparison with the data, and finally use the known 
2-loop relation between the pole mass and the MS massEJ to convert to the 
MS mass rnbijrib)- 

4 Results and discussion 

There are restrictions on the value of n that can be chosen for the moments 
(|l|) . Comparison of the NNLO threshold approximation of the cross section up 
to order a^ with the 'exact' resultLj shows that the threshold approximation 
provides a reasonable approximation as long as one is less than 2 GeV away 
from threshold. This restricts n > 6, conservatively. Requiring that all scales 
of the problem are larger than Aqcd restricts n < 10. 

The experimental moments are obtained from the known masses and decay 
constants of the T resonances below the open bb threshold. The continuum 
cross section is parametrized by the constant value 0.4 ± 0.2. (The asymptotic 
value for E -^ oo is 1/3.) For n = 8 the continuum cross section contributes 
only about 10% to the experimental moment integral. 

In Fig. nl we show the result for the 8th moment as a function of the quark 
mass. In the upper figure, the result is plotted as a function of the pole mass 
nib, in the lower figure as a function of the PS mass m(,^ps(2 GeV). For a given 
moment n, the integral (nj) implies the characteristic energy scale E ^ nib/n. 
The logarithms that appear in the PNRQCD integrals then suggest that the 
scale fi is chosen as Irribj \pn. This choice is shown as the solid line in the 
figure, together with variations by a factor 1/2, 2/3, 3/2, 2. In case of the pole 



Figure 1: The 8th moment (in GeV ^^) as function of the bottom quark pole mass (in 
GeV) (upper figure) and the bottom quark PS mass (lower figure) a,t fif = 2 GeV. The solid 
curve is for the scale /i = 2mi,/Vs, the outer (dash-dotted) curves show the result, when the 
scale is varied by a factor of 2 in both directions. The inner (dotted) curves correspond to 
2/i/3 and 3/i/2. The experimental moment is given by the grey bar. 
17 






nib 


rnb{mb) 


Remarks 


This work 


4.97±0.17 


4.26 ±0.12 


Sum rules 


ppgsEi 


4.80 ±0.06 


- 


Sum rules 


H980 


4.88 ±0.10 


4.25 ±0.09 


Sum rules 


MYggli 


- 


4.20 ±0.10 


Sum rules 


JP97/98lili 


4.60 ±0.02 


4.19 ±0.06 


Sum rules, no resummation 


PY98i 




4 AA+0.Q3 

4-44_o.04 


T(1S) mass 


Msgsii 


- 


4.41 ±0.11 


lattice HQET 



Table 1: Bottom quark mass values (in GeV) obtained from NNLO calculations. 



mass as an input parameter, the scale dependence is very large. Including a 
(small in comparison) error due to the variation of as{Mz) = 0.118 ±0.003, we 
determine mb = (4.97 ± 0.17) GeV, where the scale error is estimated from a 
variation between 2 GeV and 2/i. The scale dependence is somewhat reduced, 
but still large, with the PS mass as input. Our result is 

TOh,ps(2 GeV) = (4.60 ± 0.125(scale) ± 0.03(a^) ± 0.03(exp.)) GeV, (14) 

which translates into 



mbirfib) = (4.26 ± 0.12) GeV 



(15) 



for the MS mass. We have included an estimate of (40 ± 40) MeV for the 
unknown 3- and 4-loop terms in the relation between the PS and the MS 
mass, which are needed here, because the sum rule determines the PS mass 
with a parametric accuracy of order niba'^ . The central value is stable against 
variations of the order of the moments. The large remaining scale dependence 
can be traced to the scale dependence of the theoretical prediction of the first 
bound state residue. This is discussed in more detail in Ref. El. 

In Tab. we compare our result to other bottom quark mass results, 
choosing only NNLO calculations (although NNLO may not imply a NNLO 
resummation in some cases). The first four entries all refer to the sum rule 
method with NNLO resummation. There are several differences between the 
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present and previous implementations of the NNLO nionients_.Ihc^miiniation 
of NLO logarithms of v, which has not been done in Refs.E30E3, turns out 
to be^jSmall effect; leaving it out pspuld increase our result by up tLLlOMeV. 
Refs.EjO use the 'old' value of 02^3, which is most likely incorrectEI. Again, 
the corrected value for 02 shifts the extracted quark mass downwards by only 15 
MeV. More significant differences arise in the treatment of the short-distance 
factor and the bound state pole (5-functions. For a detailed discussion we refer 
to Ref.l. 

The sum rule calculation of Refs. Ej'EZI does not include a non-relativistic 
resummation. It leaves out, in particular, the contribution from the Coulomb 
poles. This seems to be the main reason why the battom masses of JP97/98 
come out small. The large moments used in Refs. Ea^Ej are completely domi- 
nated by the Coulomb pole contribution to (|lj). For the 8th moment, we find 
that our pole quark mass would decrease by roughly 400 MeV, if we left out 
the Coulomb pole contribution. i— . 

The last two entries in the table come from other methods. Ref.c3 uses the 
NNLO result for the Coulomb pole position and compares it with the T(1S) 
mass. The error quoted in Ref. IHa is probably underestimated, because the 
scale /i^ is allowed to vary only by ±25% rather than a factor of 2. The origiri 
of the large value of the MS quark mass in Ref. c3 will be explained in Ref. u. 
The last entry refers to the first complete NNLO quark mass determination 
from the lattice, based on the B meson mass and the lattice measurement of 
the static energy of a heavy quark. 
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